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Abstract 

A number of authors have proposed stochastic versions of the Schrodinger equation, 

either as effective evolution equations for open quantum systems or as alternative theories 

^ ^ with an intrinsic collapse mechanism. We discuss here two directions for generalization of 
> 

. these equations. First, we study a general class of norm preserving stochastic evolution 

m 
O 

' equations, and show that even after making several specializations, there is an infinity of 

o ■ 

possible stochastic Schrodinger equations for which state vector collapse is provable. Second, 
PL|| we explore the problem of formulating a relativistic stochastic Schrodinger equation, using a 
manifestly covariant equation for a quantum field system based on the interaction picture of 
Tomonaga and Schwinger. The stochastic noise term in this equation can couple to any local 



O 
O 
(N 

Ph. 
< 

(N 



^ ■ scalar density that commutes with the interaction energy density, and leads to collapse onto 



spatially localized eigenstates. However, as found in a similar model by Pearle, the equation 
predicts an infinite rate of energy nonconservation proportional to S^{0), arising from the 
local double commutator in the drift term. 
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I. Introduction 

The measurement problem is widely perceived as the greatest difficulty in the interpre- 
tation of quantum mechanics: how, without invoking a separate realm of classical measuring 
devices, can one rule out superpositions of macroscopically distinct states, as in the famous 
Schrodinger's cat paradox? To answer this question, a number of authors have suggested 
modifying the usual Schrodinger equation so as to eliminate such superpositions at large 
length scales, while retaining the standard quantum results for microscopic systems. The 
result is a modified Schrodinger equation containing extra terms, including stochastic terms 
which reproduce the probabilities of measurements [1-7] . 

In a parallel development, other researchers have derived effective equations to describe 
systems evolving in contact with an external environment. These effective equations also 
take the form of stochastic Schrodinger equations, of a form very similar to those posited in 
response to the measurement problem [8-14]. 

One example of such a modified equation is the quantum state diffusion (QSD) equation 
of Gisin and Percival [12], which has the form 

i#) = -iHmt+Y. - \LiLk - mt+Y,{Lk-{U)mdik . (i) 

Here the Lindblad operators Lk [15] represent the effects of the environment, H is the Hamil- 
tonian, and the stochastic differentials represent independent complex Wiener processes 
with vanishing ensemble averages or means (i.e., M[d^k] = 0), that obey the Ito stochastic 
calculus 

dCjdik^ dt5jk , d^jd^k^dtUjk , dtd^k ^ ■ (2) 

Equations (1) and (2) define an Ito stochastic differential equation; in manipulations using 
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the Ito differential d, one must use tlie modified cfiain rule d{AB) — dAB + AdB + dAdB. 
[In the "standard form" of the QSD equation given by Gisin and Percival, the symmetric 
complex matrix Ujk is zero. However, Wiseman and Diosi have recently shown that the 
most general Ito stochastic unraveling of the Lindblad evolution has ujk nonzero, with the 
matrix norm \\u\ \ bounded by unity [16]. Hence we will keep Ujk nonzero in setting up the 
general framework for our discussion, only dropping it later on.] While the dynamics of 1-0) 
can be extremely complex, there is a tendency for the state to localize onto eigenstates of 
the Lindblad operators L^. Of course, the competing influences of different L^, or of the 
Hamiltonian H, can prevent this localization from taking place. Note also that Eq. (1) is 
nonlinear in this will in general be necessary for such an equation to preserve the norm 
of the state. 

While we have presented this as an effective equation, arising due to the effects of an 
external environment, one can postulate an exactly similar equation in which the noise 
is considered fundamental. Percival has proposed such an equation with localization onto 
energy eigenstates, which he calls Primary sta,te diffusion (or PSD) [17]. Other such equations 
have been proposed by Pearle, by Ghirardi, Rimini and Weber, by Diosi, by Ghirardi, Pearle, 
and Rimini, and by Hughston [1,3,5,7]. A survey of their properties has recently been given 
by Adler and Horwitz [18], who give a detailed discussion of the conditions for the dynamics 
of Eq. (1) to lead to state vector collapse. 

Our aim in this paper is twofold. First, we examine the extent to which a stochastic 
dynamics such as Eq. (1) can be kept in its most general form, subject to the requirement 
that it should still lead to state vector collapse. This forms the subject matter of Sec. II, 
where we show that there is an infinite parameter family of stochastic equations for which 
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state vector collapse is provable. 

Our second aim is to explore the well-known problem that all equations with the structure 
of Eq. (1) are nonrelativistic. They are designed to mimic measurement, and they almost 
all contain a distinguished frame which takes the role of the rest frame of the measuring 
device. Since in standard QM measurements take effect instantaneously on the state vector 
of the entire system — ultimately, on the entire universe — it has been very difficult to find a 
covariant theory of measurement. In Sec. Ill we study a local generalization of Eq. (1) which 
can be written in manifestly covariant form, based on the "many fingered time" Tomonaga- 
Schwinger generalization of the Schrodinger equation. (For previous related approaches to 
this problem, see e.g. [19-23].) The generahzed equation, like its nonrelativistic counterparts, 
causes the values of certain quantities (such as the center of mass of a measuring meter) 
to localize. However, there are difficulties with energy conservation arising from the local 
structure of the stochastic terms. 

II. Generalized Stochastic Equations 

II. 1 General Framework 

We begin by giving a general framework for the basic QSD equation of Eq. (1). Consider 
the stochastic state evolution 

\d^P) = a\^P)dt + Y,^Pk\i^)dik , (3) 

k 

with d^fc independent complex Wiener processes as in Eqs. (1) and (2), and with a and (5^ 
the operator coefficients of the drift and stochastic terms respectively, which can also have an 
exphcit dependence on the state \'^). The condition for the norm of the state to be preserved 
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= d{ij\i}) = {diP\iP) + {iP\diP) + {diP\diP) . (4) 

Substituting Eq. (3) and its adjoint, and using Eq. (2) to simplify the quadratic terms in 
the Ito differentials, this becomes 

= dt{^\a + dt + + J2^daml\^) + dU^lW)] ■ (5) 

k k 

Since d^k and d^ are independent, Eq. (5) requires that the coefficients of dt, d^k, and d^ 
vanish independently, giving the conditions 



(6) 



={ij\(3k\ij) , all k . 
Letting be a set of general (not necessarily self-adjoint) operators, and H — W and 
i^'^ be arbitrary self-adjoint operators, the general solution to the conditions of Eq. (6) 
takes the form 



Pk - {ip\Lk\ip) , 

(7) 

> H; Hi. . 

k 

with the operators K and Lk still allowed to have an explicit dependence on the state vector 
It is convenient for what follows to introduce the definitions 

(6) =(V'|o|V') , 

(8) 

AO =6 - (O) , 

where O is an arbitrary operator. Then Eq. (7) can be written in somewhat more compact 
form as 

Pk=^Lk , 

(9) 



a , 

2 

k 



Equations (3) and (7)-(9) give the general form of a norm-preserving stochastic extension of 
the Schrodinger equation. Equation (1) clearly has this general form, with the specific choice 
K — I X]fc((-t'fc)-£'jk — L\{Lk))^ for which {K) = so AX = K. Usually, in applications of 
the QSD equation it is assumed that the Lindblads have no dependence on the state 1-0), 
but we will find it useful to keep open the possibility that they do have a nontrivial state 
dependence. 

To analyze convergence properties implied by this equation, we shall need formulas for 
the evolution of the expectation (O) and the variance V[0\ = ((AO)^) = (O^) — (O) of 
a general operator O. Using Eq. (2) and the Ito extension of the chain rule, together with 
Eq. (3) and its adjoint, and (in the calculation of dV) imposing the normalization constraints 
of Eq. (9), we find after some algebra the results 

d{d) ={^ + a^d + da + Y^ PlOPk)dt 

k 

+Y;^[dc,{d$,)+dcMo)] , 

k 

«fl6)=[(A6f + ^a6> (10) 

+{at(A6)^ + (AOj^d + J^/3t(A6)2A> 

k 

-2 0i^d){AdPk) - 2Re( J] {AdPk){/^dPi)uki)]dt 

k kl 

+j2m{Adr^)+dcMi^on ■ 

k 

In applying Eq. (10), we shall have occasion to take its mean over the Ito process. Since the 
stochastic expectation or Ito process mean M[ ] obeys 

M[dCkS] = M[dCkS] =0 , all A; , (11) 

for a general Hilbert space scalar S, the terms in Eq. (10) involving d^k and d^^ drop out in 



the mean, giving 



M[d{d)]=M[{^+a''6 + da + J2PkOPk)]dt , 



dt 

k 



M[dV[d]] =M[(AO^ + ^AO) (12) 
at at 

+(at(A0)2 + ^AOfa + J2Pli^0f^) - 2 (/^l AO) (AO/3,) 

k k 

-2Re{J2{^OPk){AdPi)uki)]dt . 

kl 

Clearly, these equation take the same form if O is replaced everywhere by any function 
F[0], since this simply defines a new operator F that replaces the dummy operator O. In 
the next two sections we shall argue that for the evolution given by Eq. (3) to converge to 
an eigenstate of O, we must have M[(il^[0]] < , with equality only for ((AO)^) = 0, and 
shall demonstrate this for a particular special class of equations. 

II. 2 Specialization 

We shall now introduce some simplifying specializations, which as we shall see, still leave 
an infinite parameter class of stochastic Schrodinger equations, for which state vector re- 
duction to eigenstates of the operator O is provable. First of all, let us restrict ourselves to 
the case in which O is a self-adjoint observable, which we assume to have no explicit time 
dependence, so that = 0^ , dO /dt = 0. Secondly, we now take the complex matrix ujk 
of Eq. (2) to be zero, and we speciahze the choice of the operators a and which satisfy 
the normalization constraints of Eq. (9), as follows: 

(i) We take the operator K to be zero, so that the constraint of Eq. (6) is satisfied as an 



7 



operator relation 

a + «t + ^/3t4 = , (13) 

k 

which as in Eq. (9) imphes that 

k 

(ii) We take H to be an operator that commutes with O, and take all of the Lf. to be 
functions solely of the operator O, so that they also commute with O, 

[H,d]^0 , Lk = Lk[d]^[Lk,d]^0 . (15) 

Together with Eqs. (7) and (13), these specializations imply that a and Pk all commute 
with O, as well as with any function F[0] solely of the operator O, 

(16) 

[a,F[d]]=Q , 0k,F[d]]=O , 
With these specializations, Eqs. (12) for the time derivatives of the stochastic mean of 

the quantum expectation of a function ^[O], and of the stochastic mean of the variance of 

O, simplify dramatically. Since a and (3^ commute with O, as well as with any function 

F[0], we have 

k k 
{a^F[d] + F[d]a + PkFWk) =(F[0] [a + a^ + J^ ^kPk]) = , 

k k 

where we have used the operator constraint of Eq. (13). Also, since (AO) = 0, we have 

(a64) = (a6(4-(4))) = (a6l,) , (18) 
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and when O is self- adjoint, we have {fi\.lS.O) — {ISD^^^ . Thus, what remains of Eq. (12) is 

M[rf(F[O])]=0 , 

^ ^ (19) 

M{dV{0'^=-2M'^\{t\0L^{0\)\^\dt , 

k 

with F\0\ any function solely of the operator O. 

II. 3 State Vector Reduction 

We shall now show that the stochastic dynamics, as specialized in the preceding sub- 
section, implies state vector reduction to eigenstates of O (assumed nondegenerate), with 
probabilities given by the Born rule in terms of the initial wave function. We shall need one 
further assumption beyond those introduced above, namely, that the scalar valued function 
/ of O defined by 

/[d]^5]|(A6L,[d])|^ (20) 

vanishes if and only if ((AO)^) vanishes. One simple way to achieve this is to take to 
have the form 

TV 
n=0 

with > for at least one value of k; note that here we are using the freedom, remarked 
on above, to allow the Lindblads to have an exphcit dependence on the state vector. This 
implies that for this value of k, 

N 

(AOL,[0]) = 5^4"\[(AO)2]"+^)>cf((AO)2) , (216) 

n=0 

and so the vanishing of f[0] imphes the vanishing of ((AO)^). This still leaves an infinite 
parameter freedom in the construction of the Lfc. A second specific example of a f[0] with 
the needed property is given in Sec. III. 3 below. 



A general condition for f[0] to have the needed property can be formulated by rewriting 
Eq. (20) as 

f[d] ^ {^\AdP[d]Ad\^) , P[0]^^4p]|V')(V'|^lP] , (22) 

k 

with P[0] by construction a positive semidefinite operator. If the Lindblads were all 
unity, P would be proportional to the projector \'il>){il>\, and since is orthogonal to the 
state AO|V'), one would have f[0] = 0. In order for f[0] to have the needed property, it is 
necessary for the Lindblads to introduce enough distortion of the projector for P[0] 

to make a strictly positive contribution to Eq. (22), in which case the vanishing of f[0] 
requires the vanishing of the state AO\ip), or equivalently, the vanishing of ((AO)^). This 
formulation of the condition on f[0] suggests that in the generic case, it is natural for it to 
have the needed property. 

We can now proceed with a convergence proof, following the presentation given by Adler 
and Horwitz [18] (see also [1,4,7]). Integrating the second hue of Eq. (19) with respect to t, 
we get 

M[V[d]]{t) = M[V[d]]{0)-2 [ M[f[d]]{t)dt . (23) 

Jo 

Since both V and / are nonnegative, Eq. (23) implies that the integrand M[/[0]](t) must 
vanish as t — > oo, since otherwise the right hand side of Eq. (23) would become negative 
at large times. This in turn implies that f[0]{t) vanishes as i — > oo except on a set of 
probability measure zero, which by the assumption introduced following Eq. (20) implies that 
the variance V[0] — ((AO)^) vanishes as i — > oo except on a set of probabihty measure zero. 
Thus, when O is nondegenerate, the state vector reduces to a pure state. Now integrating 
the first line of Eq. (19) with respect to t, taking the function F[0] to be a projector 11^ on 
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the ith eigenstate of O, we get 

M[(n,(oo))] = M[(n,(o))] = (n,(o)) . (24) 

The left hand side of Eq. (24) is just the probabihty that the stochastic process settles at 
i = oo on the ith eigenstate of O, while the right hand side of Eq. (24) is the probability 
amplitude squared for the £th eigenstate to occur in the initial state vector Thus, as first 
proposed by Pearle [1], the first line of Eq. (19) — which states that the stochastic process for 
{F[0]) is a martingale — imphes that the state vector reduction implied by the second fine 
of Eq. (19) obeys the Born probability rule. 

To compare what we have done to the analyses of Hughston and of Adler and Horwitz 
[7,18], those authors consider the energy-driven case in which the operator O = H, and in 
which one (at least) of the AL^ is simply taken as AH, corresponding to the case where 
the sum in Eq. (21) consists only of the n = term. Note that in this case, it makes no 
difference whether we take — H or we take = AH, since either gives AL^ — AH. 
When n > terms are present in L^., this distinction is important, and plays a role in our 
example showing that there are more general stochastic equations that still allow one to 
prove state vector reduction. 

III. Relativistic Stochastic Equations 

Because Eqs. (1) and (3) involve a universal time variable t at all spatial points, they are 
clearly nonrelativistic. This is evident from Eq. (2), which states that the same Ito stochastic 
differential is present everywhere in space, giving Wiener processes at space-like separated 
points that are totally correlated. In this section we explore the possibility of extending 
Eq. (1) into an equation with local Wiener processes, which can then be generalized to 
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manifestly covariant form. We shall work henceforth with a relativistic quantum field theory, 
rather than with a nonrelativistic quantum mechanical system, and thus will seek to modify 
the Schrodinger equation for this field system to a stochastic Schrodinger equation analogous 
to Eq. (1), in a manner that preserves relativistic covariance. For closely related work, from 
which our analysis differs in some details, see Pearle [19], Ghirardi, Grassi, and Pearle [20], 
and Diosi [23]. 

III.l The interaction picture 

Suppose we choose a particular Lorentz frame with coordinates t, x, and define a state 
vector l'^) for a field system at time t. This state evolves to a new state at time t + dt 
according to the Schrodinger equation 

\dip{t)) - -iH\ip{t))dt . (25) 

The Hamiltonian H is the integral of a Hamiltonian density H{x) over a constant time 
surface, 

J d^xH{x) . (26) 

Neither the Hamiltonian H nor the Hamiltonian density H{x) are Lorentz invariant, and 
the Hamiltonian densities at points x and y do not commute, 

[H{x), H{y)\ = -iVs5\x - y) • (P{x) + P{^) , (27) 

with P[x) the momentum density. These facts make it difficult to directly extend Eq. (25) 
into a stochastic Schrodinger equation in a manner consistent with Lorentz invariance. 

As a first step in avoiding these problems, let us switch to the interaction picture. (Our 
use of this is heuristic and ignores mathematical issues of the existence of the interaction 

12 



picture, as discussed e.g. in [24].) We write the Hamiltonian density H{x) as a sum 

H{x)^Ho{x) + H,^,{x) , (28) 

where Hq{x) is the free-field Hamiltonian density and Hint{x) is the interaction Hamilto- 
nian density. Unlike H{x) as a whole, Hint{x) is a relativistic invariant in theories without 
derivative couplings, and it commutes with itself at different points, 

[^int(f),^int(y)]=0 . (29) 

Let U be the unitary time evolution operator for the free-field Hamiltonian, 

U = exp{iHot} . (30) 

If lipsit)) is the state at time t in the Schrodinger picture, the state in the interaction picture 
is \ip{t)) = U\ips{t))- We similarly replace the Schrodinger picture field operators (e.g., 0, tt) 
with interaction picture operators (e.g., (/)(t) = U(l)U\n{t) — UnW). If we express the 
interaction Hamiltonian density Hiat{x), which is a function of the field operators at the 
point X, as a function of the interaction picture field operators, the state then obeys the 
simple evolution equation 

\dij{t)) ^ -iHi^tmt))dt , (31) 

where 

^int = j d^xHi^tix) . (32) 

One must now remember that operators that were time independent in the Schrodinger 
picture acquire a time dependence governed by the free Hamiltonian Hq. 
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So far this discussion has been restricted to constant-time surfaces in a single Lorentz 
frame, in which a fixed time step dt is taken simultaneously at all spatial points x, and 
so Eq. (31) is still not Lorentz invariant. We now follow Tomonaga and Schwinger [25,26] 
(see also Matthews, KroU, and Dyson [27-29]) in generalizing Eq. (32) into a /oca/ evolution 
equation. Consider a spacelike surface a with local coordinates x, on which the state of 
the underlying quantum fields is described by a Fock space state vector \ip{a)). Instead of 
advancing the whole spacelike surface a, we instead move the surface forward (i.e., in the 
normal direction) by an increment dt{x) only in the vicinity of a single point x, distorting 
the surface cr to a new spacelike surface a'. Under this evolution, the state vector \ip) evolves 
to a state vector l'^) + rf^jV'), with the change in the state vector given by 

dg\^lj) ^ -iH,nt{xmdt{x) . (33) 

The change in the state vector resulting from advancing the entire surface is then 

|#) = J d^Sds\iP) . (34) 

Since the Hij^t{x) at all points commute, the order in which the spacelike surface is advanced 
is immaterial and so the right hand side of Eq. (34) can be unambiguously integrated, 
which for constant-time surfaces with dt{x) = dt recovers the original interaction picture 
Schrodinger equation of Eqs. (31)- (32). From the viewpoint of constructing a stochastic 
generalization, the local form of the interaction picture evolution equation given in Eq. (33) 
has three advantages: It is readily put in manifestly covariant form, it involves only the 
Lorentz scalar operator density ^int(a;), and this operator commutes with itself (and with 
other easily constructed scalar densities) at spacelike separations. 

III. 2 The local norm-preserving stocheistic equation 
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Let us now replace the local unitary evolution equation of Eq. (33) with a new equation 

d^lV') = a{x)\ilj)dt{x) + P{x)\^|J)dC{x) , (35) 

in which we take the coefficient functions a{x) , P{x) and a{y) , P{y) to mutually commute 
for all X , y, so that no noncommutativity problems are encountered when we compound 
evolutions for different values of x. Here d^{x) is a complex stochastic differential variable 
defined at each point x, which has zero stochastic mean (i.e., M[d^{x)] — 0), and which 
obeys the local Ito calculus 

dCix)day) = S\x-y)dt{x) , di{x)di{y) = dt{x)di{y)=Q . (36) 

The spatially integrated form corresponding to Eq. (35) is 

In analogy with our discussion of Sec. II. 1, we can now determine the conditions on the 
coefficient functions a{x) and (5{x) for Eq. (37) to preserve the norm of the state, 

d^x (-01 [di{x) + a{x)'^ + /3(f )^/3(f )] |'0)cit(f ) (3g) 

+ j d''xmm^i^)dc{s) + {i^\mwdm] ■ 

Since d^*{x), d^{x), and dt{x) are linearly independent, the normalization of the state 
(-01-0) — 1 is preserved if and only if for all x we impose the conditions 

0={ij\[a{x) + a{xy + p{xyp{x)]\ij) , 

(39) 

O^i^l^lPixM . 

Evidently, if we were to replace a in Sec. Ill by d^, then Eq. (39) could be viewed as a 
local version of Eq. (6), with x playing the role of the index k. Imposing the normalization 
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conditions, and specializing henceforth to dt{x) = dt and flat spacehke surfaces a, we flnd 
the following local version of Eq. (10), 



d{0) ={dO) + J d^x[{a\x)0 + Oa{x) + (3\x)0(3{x))]dt 

+ [ d'x[dax){dm) + dC{x)0^{x)d)] , 



dV[d] ^{AOdO + dOAO) 

+ J d^x[{a\x){Adf + {Adfa{x) + p\x){Adfp{x)) - 20\x)Ad) {AO p{x))]dt 

+ J d'xmx){{Adfm)+de{s)0\x){Adf)] . 

(40) 

Instead of working with the most general form of the normalization condition, we shall 
speciahze (as we did in Sec. II. 2) and satisfy Eq. (39) by taking a{x) and P{x) to have the 
form 

a{x) = - iHint{x) - ^/3(f)^/3(f) , 

2 (41) 

m=^S{x) , 

with S{x) any local Lorentz scalar operator that commutes with Hiai{x). We shall further 
assume S{x) to be self-adjoint. Additionally, we shall assume that the operator O is self- 
adjoint and has no intrinsic time dependence in the Schrodinger picture, so that in the 
interaction picture its time dependence is given by 



f = .[i/..6l . 
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(42) 



With these speciahzations, Eq. (40) can be rewritten after a httle algebra as 



d{d) 



m,d]-^ j d'x[s{x),[s{x)A]) 



dt 

+ j d''x[di{x){0/\S{x)) +dC{x){^S{x)d)] , 

(43) 



dv[d] = 



dt 



{i[H,{Adr]-^ j d'x[S{x),[S{x),{Adn)-2 J d'x\{AdAS{x))\ 
+ J d^x[d^{x){{AdyAS{x)) + dC{x){AS{x){Adf)] . 

Although not needed for our purposes, by using the fact that d^xdt and S{x) are Lorentz 
scalars, the stochastic and drift terms in Eq. (43) can be readily written in manifestly 
covariant form. The corresponding covariant transcription of the Hamiltonian evolution 
terms is given in Matthews [27] and KroU [28] . 

III. 3 Reduction for local density eigenstates 

Let us now apply the above formulas to discuss state vector reduction to local density 
eigenstates, giving a relativistic generalization of the localization models discussed in [3-5]. 
Let us make the specific choice 

S{x) ^ CH,^t{x) , (44) 

which obviously satisfies the commutativity conditions 

[a{x), am = [aix),m] = [/9(f), ^(y)] = [a(f), ^int(^] = [m, HUv)] = , (45) 

for all spacehke separated points x , y. In field theories like the Standard Model, in which 
all mass comes from spontaneous symmetry breaking, the mass terms arise from ifint(f), 
and so for bulk matter we are effectively taking S to be the local mass density operator, 
multiphed by a scale factor C. 
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As a concrete illustration of how Eq. (43) can lead to state vector reduction and local- 
ization, let us consider the simplified case of an apparatus connected to a pointer with two 

— * — * 

macroscopic states specified by two values Xi , X2 of the pointer center of mass variable 

X = -^P""^*^" (46) 

r . , d^xs(x) 

J pointer ^ 
A — * 

We shall apply Eq. (43) to this system, taking O — X. By Eq. (45), the double commutators 
[S{x), [S{x), O]] and [S{x), [S{x), (AO)^]] both vanish, but in general the commutators [H, O] 
and [H, (AO)^] are nonzero. However, if we take the two macroscopic pointer positions to 
be degenerate in energy, then the commutators involving H vanish within the degenerate 
two-state subspace. Taking the stochastic mean M[ ] of Eqs. (43), we then find within the 
two-state subspace the simplified equations 

M[d{X)] =0 , 

r ^ (47) 

M[dV[X]\ = - 2 / d^xM[\{AXAS{x))\^]dt . 

We have here exactly the same structure as we found in Eq. (19) above, and the function 
f{X) = J d^f|(AXA^(f))p is easily seen [c.f. the final line in Eq. (48) below] to obey the 
condition that the vanishing of f{X) implies the vanishing of ((AX)^). Hence the same 
argument as was used in Eqs. (23) and (24) proves that an initial superposition of the two 
center of mass eigenstates reduces to either the state with X = Xi or the state with X = X2, 
with respective probabilities given by the amplitude squared to find the initial state in the 

— * 

respective X eigenstate. 

From Eq. (47), we can estimate the reduction rate T as follows. Writing \ip) = \Xi) cos 9+ 
1X2) sm9, and assuming that the states \Xi) , 1X2) differ sufficiently for us to approximate 
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(48) 



that {Xi\S{x)\X2) — 0, we have after a short calculation 

V[J(] ={{AJlf) = sin^ ^cos' e{Xi - X^f , 
{AXAS{x)) =sin2^cos2^(Xi -X2)((Xi|5(f)|Xi) - {X2\S{x)\X2)) , 
\{AJIAS{x))\'' ^sm^9cos''e{X,-X2f\{Xr\S{x)\Xi) - {X2\S{x)\X2)\'' 

{x,-x2y 

Thus Eq. (47) becomes 

^^'^^"'''Ir""'^^^ " d'xM[sm''ecos^e]\{Xi\S{x)\X^) - {X2\S{x)\X2)\'' , (49) 

from which we see that, up to numerical factors of order unity, the reduction rate is given 

by 



r~ /" d'x\{x,\six)\Xi) - {X2\s{x)\X2 

(50) 

~CM Mass Density]^ . 

J pointer 

For a pointer containing ~ 10^^ nucleons of mass M ~ 1 GeV and volume V ~ 10~^^ cm^, 
the estimate of Eq. (50) becomes 

r ~ C^NM^V'^ , (51) 

which gives a reduction rate T > 10^ sec~^ (corresponding to a collapse time faster than 
characteristic observational time scales) for C > (10^ GeV)~^. This corresponds to a mass 
scale at roughly the geometric mean between the Planck mass and a nucleon mass. Thus, 
in contrast to the energy driven model [7,12,18] for state vector reduction, where the mass 
scale for the coefficient of the noise terms is Planckian, in the local version discussed here 
the mass scale for the noise terms is much below the Planck scale, but still large compared 
to elementary particle masses. 
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III. 4 Energy nonconservation 

Except for the special case of stochastic equations in which the Lindblads are taken to be 
operators that commute with the Hamiltonian (including the Hamiltonian itself), stochastic 
modifications of the Schrodinger equation lead to energy nonconservation, as has been noted 
in the papers of Ghirardi, Rimini, and Weber, and of Pearie [1,3,19]. Let us examine this 
issue in the context of the relativistic model discussed above. For any operator O, the 
stochastic expectation M[ ] of the first formula in Eq. (43) is 

1 



M[d{0)] = M 



mA-^ I d'S[S{x),[S{x)M]) 



dt , 



(52) 



which when applied to the Hamiltonian (i.e., taking O — H) gives for the mean rate of 
energy nonconservation 



M 



'd{H)' 




dt 





J3-? 



(53) 



In typical field theory models, the double commutator appearing in Eq. (53) is not only 
nonzero, but as first noted by Pearie [19] is proportional to 5^{Q) and thus is infinite. For 
example, taking a Dirac field model with 

j d''xi^\x)[i-^a ■ V + I34){x)]i^{x) ^Ho + Hint , (54) 

with 0(x) an external scalar field with nonzero vacuum expectation, and choosing 

S{x)=CHint{x) = Cij\x)f3<P{x)'iP{x) , (55) 

one has 

[S{x), [S{x), H]] = S^{0)C^<|)''{x)i-^[^p''{x)d ■ Vt^ix) - ^^p\x) ■ a^Pix)] . (56) 
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Similar results are found in scalar meson field theory models, and appear to be generic. 
Moreover, except for special choices of S{x) (see, e.g. [19-21]), the coefficient of 5^(0) is a 
nontrivial operator and not a constant. The S^{0) singularity is a direct result of the local 
derivative structure of the drift term, and we have not found a mechanism to cancel it within 
the standard stochastic differential equation and quantum field theory framework discussed 
here. 

IV. Conclusions 

We have presented two generalizations of stochastic Schrodinger equations for state vec- 
tor collapse. First, we have shown that there is an infinite parameter family of such equations 
for which one can prove state vector collapse with probabilities given by the Born rule. Sec- 
ond, we have given a relativistic stochastic equation which can be made manifestly covariant, 
and which produces localization onto mass density eigenstates. This produces spatial local- 
ization for superpositions of macroscopically distinct system states; to give rapid enough state 
vector localization in plausible experimental setups, the scale mass governing the stochastic 
terms must be considerably smaller than the Planck mass. The local equation has the defect 
that it leads to a divergent rate of energy nonconservation in generic field theory models, 
indicating that new ideas will be needed to achieve a satisfactory relativistic state vector 
collapse model. 
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